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ABSTRACT 


Several  software  enhancements  to  the  Towed  Array 
Performance  Prediction  System  (TAPPS)  are  described.  A 
technical  discussion  of  the  changes  is  given,  and  the 
implementation  and  testing  of  the  software  modifications  on 
the  Nova  820  computer  are  documented.  The  latest  software 
is  denoted  as  TAPPS  Version  1.2. 

An  interim  improved  treatment  of  bottom  loss, 
denoted  as  interim  Bottom-Loss  Upgrade  (BLUG)  is  discussed. 
-A  set  of  six  new  bottom-loss  classes  and  associated  bottom- 
loss  curves  for  the  Indian  Ocean,  and  new  bottom-loss  curves 
in  the  other  ocean  areas  are  presented.  These  changes  have 
been  fully  integrated  with  the  TAPPS  propagation-loss  model 
(FACT),  as  well  as  the  TAPPS  Noise  Model. 

A new  surface-duct  propagation-loss  model  is  a 
refinement  and  extension  of  Labianca's  virtual  mode  results. 
The  results  of  Spofford  and  Ryan  had  previously  been  used  to 
develop  a replacement  model  for  the  Clay  Surface-Duct  Mod^l 
within  FACT.  The  implementation  of  these  results  in  TAPPS 
are  discussed . ^ 

\ 

The  TAPPS  shipping-noise  data  base  has  been 
extended  to  include  the  Indian  Ocean  and  the  capability  to 
handle  data  in  the  southern  hemisphere  has  been  provided 
for. 


A Rough  Surface-Loss  Model  in  FACT  and  other 
miscellaneous  TAPPS  software  modifications  are  also 
discussed. 
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Section  1 
INTRODUCTION 

The  Towed  Array  Performance  Prediction  System 
(TAPPS)  was  developed  under  the  Naval  Ocean  Research  and 
Development  Activity  (NORDA)  sponsorship,  during  the  years 
1976  to  1979.  It  was  designed  as  a seagoing  minicomputer 
system  to  provide  Fleet  personnel  with  estimates  of  towed 
array  performance  and  to  advise  tactical  decision  makers 
concerning  the  selection  of  such  parameters  as  ship  speed, 
course,  array  depth,  and  search  frequency. 

TAPPS  was  developed  using  state-of-the-art 
modeling  techniques  and  data  bases,  as  they  existed  in  1977. 
However,  during  the  last  four  years  new  acoustical  models 
and  oceanographic  data  bases  were  developed  by  the  scien- 
tific community.  Specifically,  an  interim  Bottom-Loss 
Upgrade  (BLUG)  and  a new  Surface-Duct  Propagation-Loss  Model 
(SDUCT)  were  developed  and  implemented  at  the  Fleet  Numer- 
ical Oceanography  Center  (FNOC)  in  Monterey,  California. 
These  enhancements  became  an  integral  part  of  the  CDC  6600 
version  of  the  FACT  Propagation-Loss  Model.  Since  TAPPS  also 
uses  FACT  (modified  to  run  on  the  Nova  800  computer),  the 
decision  to  include  these  enhancements  in  TAPPS  was  made. 

During  the  course  of  implementing  interim  BLUG 
and  SDUCT  within  TAPPS  several  other  areas  for  model  im- 
provement were  identified  and  included  in  the  current 
effort.  This  included  Ambient  Noise  Data  Base  Improve- 
ments, Rough-Surface  Losses  within  FACT,  and  four  program 
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"bugs".  TAPPS  with  all  ol  these  improvements  has  been 
designated  TAPPS  Version  1.2.  A previous  version  of  TAPPS 
(Version  1.1)  was  fully  documented  in  References  1 through 
7.  This  report  discusses  only  those  features  unique  to 
TAPPS  Version  1.2.  Program  listings  for  TAPPS  Version  1.2 
are  available  in  References  8 and  9. 


Section  2 
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SOFTWARE  ENHANCEMENTS  - TECHNICAL  DISCUSSION 


This  section  discusses  the  technical  issues 
related  to  the  TAPPS  software  enhancements. 

. i 

' 1 2.1  INTERIM  BOTTOM-LOSS  UPGRADE  (BLUG) 


■tlj 


I 

r • 

as 


TAPPS  Version  1.1  used  two  sets  of  bottom-loss 
classes  to  make  acoustic  propagation-loss  and  ambient-noise 
predictions.  For  low  frequencies  (below  1000  Hz),  estimates 
were  based  on  a set  of  three  curves  of  loss  versus  grazing 
angle  applied  over  specific  frequency  bands.  The  three 
curves  are  assigned  based  on  a data  base  of  bottom  classes 
(1,  3 or  4)  which  covers  most  of  the  northern  hemisphere  at 
1-degree  intervals.  These  classes  and  associated  bottom 
loss  curves  are  actually  a consolidation  of  five  classes 
developed  in  1970  (Reference  10) . For  high  frequencies 
(above  1000  Hz) , a set  of  ten  classes/curves  were  used  based 
on  Reference  11.  The  TAPPS  Version  1.1  bottom-class  data 
base  and  bottom-loss  curves  were  essentially  the  same  as 
those  used  at  FN0C,  Monterey,  except  that  no  interpolation' 
was  performed  between  1000  and  1500  Hz  to  smooth  out  poten- 
tial discontinuities  in  acoustic  predictions. 
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Further  improvements  to  the  approach  used  pre- 
viously by  FN0C  and  in  TAPPS  Version  1.1  were  developed  in 
1P8C  and  implemented  at  FN0C  (Reference  12).  First,  a new 
geographic  assignment  of  the  three  low  frequency  classes  was 
developed  as  a result  of  the  recognition  that  bottom  loss 
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appeared  to  correlate  with  sediment  thickness.  The  three 
classes  (still  designated  1,  3 and  4)  were  retained  as 
somewhat  representative  of  measured  data  from  poor,  moderate 
and  high-reflectivity  areas.  However,  as  more  bottom-loss 
data  yielded  to  interpretation  in  terms  of  geo-acoustic 
models,  and  comparisons  of  model  predictions  with  measured 
data  reinforced  these  interpretations,  a new  set  of  curves 
was  developed  based  upon  geo-acoustic  modeling  and  data 
interpretation  to  go  with  the  three  areas  of  different 
sediment  thickness.  The  high  frequency  classes/curves  were 
not  modified . 

Finally,  recent  work  by  Applied  Research  Labora- 
tories (ARL)  and  SAI  has  produced  a set  of  six  low-frequency 
bottom  classes  and  associated  bottom-loss  curves  (designated 
1-6)  for  the  :,id4an  Ocean  (Reference  13).  These  bottom-loss 
improvements  were  also  implemented  at  FNOC,  Monterey  during 
1980. 


TAPPS  has  subsequently  been  modified  to  include 

most  of  the  interim  BLUG  features  currently  available  at 

FNOC.  The  new  low-frequency  bottom-class  assignment  for  the  . 

Indian  Ocean  (Figures  1 and  2)  have  been  implemented.  Also,  \ 

; 

new  bottom  loss  versus  grazing  angle  functions  (BTMLOS, 
BTMHF  and  BTMLF)  were  written  to  include  the  new  1,  3,  4 
class  curves  (Atlantic,  Pacific  and  Mediterranean)  and  the 
1-6  class  curves  (Indian  Ocean).  In  order  to  distinguish 
between  the  two  sets  of  curves,  a preprocessor  maps  classes 
1,  3,  4 and  1-6  into  classes  1-3  and  4-9,  respectively.  An 
interpolation  routine  was  written  to  smooth  out  discontinu- 
ities between  the  low  and  high  frequency  predictions. 
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Figures  3 through  11  show  the  new  bottom  loss  versus  grazing 
angle  curves  for  classes  1-9. 

The  current  bottom  loss  curves  and  areas  must  be 
viewed  as  an  interim  improvement,  since  a more  substantial 
change  will  be  installed  at  FNOC  within  the  next  year. 

Function  BTMLOS  is  called  in  TAPPS  by  both  the 
FACT  propagation-loss  model  as  well  as  the  ambient-noise 
prediction  model.  Also,  constants  defining  critical  grazing 
angle  versus  bottom  class  were  also  modified  to  be  consis- 
tent with  the  new  BTMLOS  function. 

2.2  SURFACE- DUCT  MODEL  UPGRADE 

TAPPS  Version  1.1  used  the  Clay  model  within  FACT 
for  surface-ducted  propagation  (Reference  14) . While  the 
Clay  model  was  recognized  to  have  certain  deficiencies,  an 
efficient  alternative  was  not  available.  The  model  had  an 
over-simplified  depth  dependence  and  the  wrong  dependence  of 
leakage  from  the  duct  on  environmental  parameters. 

Following  a series  of  evaluations  which  identified  ' 

l 

the  short-comings  of  the  Clay  model,  SAI  proceeded  to 
develop  a substantially  improved  Surface-Duct  Model  based 
largely  on  the  work  of  Labianca  (Reference  15),  Spofford 
(Reference  16),  and  Ryan  (Reference  17). 

The  Surface-Duct  Model  is  designed  to  compute  the 
relative  incoherent  intensity  level  from  a point  source  in  a 
canonical  bilinear  duct.  The  approach  used  in  the  new  model 
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Figure  3:  Bottom  Loss  Curves 
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Figure  5:  Bottom  Loss  Curves  (Class  3) 
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Figure  7:  Bottom  Loss  Curves 


Figure  9:  Bottom  Loss  Curves  (Class  7) 
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is  a refinement  and  extension  of  Labianca’s  virtual  mode 
results  to  include  modes  near  cutoff  as  well  as  leaky 
modes. 


Rough  surface-reflection  losses  are  computed  in 
a separate  module  (see  Section  2.4)  as  a function  of  the 
ray's  surface  grazing  angle,  wave  height  and  frequency. 

The  integration  of  the  new  Surface-Duct  Model 
required  the  development  of  a preprocessor  to  transform  the 
surface-duct  portion  of  the  sound  speed  profile  into  a three 
point  (bilinear)  profile.  A routine  (SDINIT)  was  developed 
to  perform  this  function  as  well  as  error  checking. 

2.3  AMBIENT  NOISE  UPGRADE 

TAPPS  Vet  o ton  1.1  contains  average  shipping-noise 
levels  for  each  five-degree  square  in  the  northern  hemi- 
sphere from  0 to  75  degrees  North.  Although  provision  was 
made  to  include  the  Indian  Ocean  in  this  data  base,  this 
portion  of  the  file  was  left  empty  due  to  the  lack  of  data 
at  the  time  the  file  was  created  (1977).  Recently  several 
sources  of  noise  measurements  in  the  Indian  Ocean  have1 
become  available  (References  18,  19  and  20).  Hence,  the  ’ 

decision  *o  extend  the  TAPPS  shipping  noise  file  (ANFILE)  to 
include  the  Indian  Ocean  was  made. 

The  noise  measurements  were  for  16  separate  areas 
in  the  Indian  Ocean,  and  were  typically  time  and  frequency 
averaged  measurements  at  relatively  deep  deoths  (at  least 
300  feet).  This  minimised  the  possibility  oi  contamination 
by  surface- image  interference  and  other  dynamic  effects. 


The  data  base  was  derived  by  overlaying  on  a 
single  grid  the  basic  noise  measurement  levels,  shipping 
density  (low,  medium  and  high),  and  areas  with  depth  excess. 
Regions  with  no  depth  excess  were  denoted  by  the  bottom 
class  (low,  medium  or  high).  Based  on  all  of  these  factors, 
areas  with  no  noise  measurements  were  filled  in  by  interpo- 
lation/extrapolation. Values  were  determined  for  each 
five-degree  square  in  the  Indian  Ocean  from  10°S  to  30<>N, 
This  was  necessary  to  be  consistent  with  the  existing 
water-mass/bottom-class  files  in  TAPPS  which  cover  these 
same  locations. 

The  existing  ANFILE  was  expanded  to  include  the 
southern  hemisphere,  and  provision  to  access  the  southern 
hemisphere  portion  of  the  file  was  made  by  modifying  sub- 
routine SHPSR,  Although  only  a small  section  of  the  Indian 
Ocean  uses  the  southern  hemisphere,  future  modifications  to 
incluoe  other  areas  of  the  world  can  now  be  made  with 
minimal  effort.  Several  points  in  the  Mediterranean  were 
found  to  be  in  error  and  corrected. 

The  revised  ANFILE  is  shown  in  Figure  12. ^ The 
first  half  of  the  file  contains  the  northern  hemisphere  and 
the  second  half  the  southern  hemisphere. 

2.4  ROUGH-SURFACE  LOSS  UPGRADE 

A routine  to  compute  rough  surface  losses  was 
developed  based  on  Reference  21. 


The  surface  loss  is  computed  as: 
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Figure  12:  Shipping  Level  File  (ANFILE) 
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Figure  12  (Continued):  Shipping  Level  File  (ANFILE) 
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Figure  12  (Continued):  Shipping  Level  File  (ANFILE) 
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SURLOS  » (e-0)NB 


where 

NB  = number  of  surface  reflections 
G = (WH2/8)  ( 4ttF/ 5000)2  sin2  0 
WH  = wave  height  (ft) 

F = frequency  (Hz) 

8 = grazing  angle  at  surface 

The  parameter  G is  limited  to  a maximum  of 

loge(10),  to  limit  SURLOS  to  a maximum  of  10  dB  per  surface 

reflection. 

2.5  MISCELLANEOUS  MODIFICATIONS 

Several  minor  software  bugs  were  discovered  during 

the  course  of  software  testing.  Fixes  for  these  problems 
are  discussed  in  Section  3.5. 


2-22 


Section  3 

SOFTWARE  IMPLEMENTATION  NOTES 


This  section  summarizes  the  software  modifications 
required  to  implement  the  enhancements  discussed  in  Section 
2.  The  revised  program  listings,  load  maps,  etc.  are  not 
included  due  to  the  large  number  of  pages  involved  (Refer- 
ences 8 and  9).  They  are  available  at  either  Science 
Applications,  Inc.  (McLean,  VA)  or  at  Analysis  and  Tech- 
nology (North  Stonington,  CT) . 

A summary  of  the  source-code  modifications  is 
shown  in  Table  1.  Almost  all  of  the  TAPPS  overlays  are 
affected  by  the  software  modifications.  A further  breakdown 
and  discussion  of  the  program  modifications  is  contained  in 
Sections  3.1  to  3.5.  Tables  2 and  3 contain  the  latest 
TAPPS  master  disks  maps;  as  of  16  June  1981. 

3.1  INTERIM  BOTTOM-LOSS  UPGRADE  (BLUG)  SOFTWARE 
MODIFICATIONS 

3.1.1  TAPPS  BLUG  Software  Modifications 


INPUTOV;SR  Ocean  areas  1-4  are  defined  to  dis- 
tinguish between  the  Pacific  (1),  Atlantic  (2), 
Mediterranean  (3)  and  Indian  Ocean  (4).  Ocean 
areas  1-3  use  bottom  classes  1-3  and  ocean  area  4 
now  uses  classes  4-9.  A call  to  Function  IBCON  is 
made  to  unpack,  reset  and  pack  the  bottom  class 
variable  IB  based  on  ocean  area.  Unused  code  was 
removed  to  make  room  for  the  latest  changes. 
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s 
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s 

Keyt  M - Main  program  modifiea 

S - Subroutine  called  frctn  overlay  modified 
L - Loader  file  modified 
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AN : SR  The  array  STHC  was  increased  in  size  from 
3 to  9 to  handle  the  new  Indian  Ocean  classes  4-9. 
Values  for  STHC,  the  sin  of  the  critical  grazing 
angle,  were  redefined  to  be  consistent  with  the 
new  bottom-loss  curves  incorporated  in  Function 
BTMLOS.  The  new  critical  angles  are  defined  for 
classes  1 to  9 to  be  40,  30,  0,  30,  30,  30,  30,  0 
and  0 degrees,  respectively. 

MEASUR:SR  The  size  of  the  array  STHC  was  in- 
creased from  3 to  9 (see  AN: SR). 

MEASOV : SR  The  size  of  the  array  STHC  was  in- 
creased from  3 to  9 (see  AN:  SR). 

FACTL1 :SR  The  critical  grazing  angle  (THETMB)  in 
FACT  was  redefined  as  0.5  radians  for  all  bottom 
classes.  The  array  THETCR(3)  is  no  longer  used. 

S0N0: SR  The  size  of  the  array  STHC  was  increased 
from  3 to  9,  and  the  values  modified  to  be  consis- 
tent with  the  new  BTMLOS  (see  AN: SR).  The  same' 
changes  were  also  made  to  subroutine  SN0ISE.  \ 

TAPPS : SR  The  low-frequency  bottom  class  (packed 
in  IB)  is  no  longer  reset  in  subroutine  STRC  since 
1BC0N  now  performs  this  function  in  overlay 
INPUTOV. 

BTMLOS? SR  This  Function  was  completely  rewritten 
to  compute  the  revised  low-frequency  bottom-loss 
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curves  for  classes  1-9,  and  to  interpolate  between 
the  low-  and  high-frequency  curves  between  1000 
and  1500  Hz.  New  Functions  BTMLF  and  BTMHF  were 
written  t.o  perform  the  low-  and  high-frequency 
calculations,  respectively. 

I BOON : SR  A new  Function  IBCON  was  written  to 
transform  the  low-frequency  bottom  class  obtained 
from  the  TAPPS  bottom-class  data  files.  For  ocean 
areas  1,  2 and  3,  classes  1-5  are  transformed  to 
1-3;  and  for  area  4 classes  1-6  are  transformed  to 
4-9.  The  new  low-frequency  class  is  then  repacked 
with  the  original  high  frequency  class  in  IB.  The 
transformation  was  required  to  give  BTMLOS  a way 
to  distinguish  between  the  two  sets  of  bottom-loss 
curves. 

DBLIB:SR  (BNGET)  The  arrays  THCR  and  CTHC  were 
increased  in  size  from  3 to  9 to  accommodate  the 
new  bottom  classes  4-9. 

N0ISE:SR  The  arrays  THCR  and  CTHC  were  increased, 
in  size  from  3 to  9 to  accommodate  the  new  bottom  \ 
classes  4-9.  New  critical  angles  were  defined  to 
be  consistent  with  the  new  Function  BTMLOS  (see 
AN: SR) . 

ANGSCH : SR  Debug  printout  was  added  to  aid  in 
model  checkout. 
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TWDPT:SR  Unused  code  was  deleted  to  make  room  for 


the  latest  changes. 

INDABLP/ INDBBLP  The  Indian  Ocean  bottom-class 
data  files  were  updated  with  new  low-frequency 
bottom  classes  (1-6) . The  high-frequency  classes 
were  not  changed.  The  file  coverage  extends  from 
10  degrees  South  to  30  degrees  North  in  the  Indian 
Ocean . 

BLUG  Support  Software 

LSCBLP: SR  This  routine  modifies  the  Indian  Ocean 
bottom-class  data  files  INDABLP  and  INDBBLP.  It 
accepts  as  input  two  ASCII  data  files  INDAINP  and 
INDBINP,  which  contain  the  new  low-frequency 
bottom-class  geographic  assignment  by  one-degree 
squares.  Also  inputted  are  the  original  bottom 
class  files  INDABLP  and  INDBBLP,  which  contain 
packed  low-  and  high-frequency  entries.  The 
routine  outputs  new  data  files  INDAOUT  and  INDBOUT 
with  the  new  low-frequency  classes,  packed  with 
the  original  high-frequency  classes. 

LSDBLP : SR  This  routine  reads  and  displays  any 
TAPPS  bottom-class  file  (e.g.  , INDABLP). 

LSBTM ; SR  A test  driver  for  subroutine  BTMLOS. 
The  routine  lists  bottom  loss  versus  grazing  angle 
for  each  class  and  for  selected  frequencies. 
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LSDBLOS:SR  A test  driver  tor  subroutine  BTMLOS. 


The  routine  graphically  displays  bottom  loss 
versus  grazing  angle  for  each  class  s.nd  for 
selected  frequencies. 

LSLIBrSR  A utility  library  for  character  input 
and  manipulation.  It  contains  routines  MOVEIT, 
KCOMP  and  C INPUT. 

3.2  SURFACE-DUCT  MODEL  SOFTWARE  MODIFICATIONS 

3.2.1  TAPPS  Surface-Duct  Software  Modifications 


FACTL1 :SR  The  old  Surface-Duct  Model  (Clay)  equa- 
tions were  removed.  Subroutine  SDINIT  is  called 
to  see  if  their  is  a surface  duct.  IML  and  IMLP 
are  set  to  zero  if  their  is  no  surface  duct,  or  an 
error  condition  exists  which  would  prevent  the 
successful  completion  of  the  Surface-Duct  Model. 

FACTL2 :SR  The  old  Surface-Duct  Model  contribution 
and  volume  absorption  loss  were  deleted  and 
replaced  with  a call  to  a third  FACT  overlay  \ 
(FACTL3).  * 

FACTL3 :SR  The  new  FACT  Surface-Duct  Model.  It 
first  calls  SDINIT  to  compute  the  Surface-Duct 
Model  input  parameters.  Next,  it  calls  SDUCT  if 
IMLP  is  greater  than  one  (if  a duct  is  present  and 
no  error  condition  exists).  Next,  it  adds  the 
surface-duct  contribution  and  volume  absorption 
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to  the  previously  computed  TL  array  which  includes 
signal  intensity  levels  from  all  other  propagation 
paths  other  than  surface  duct.  The  TL  array  is 
converted  to  dB  and  then  to  integers  to  save  space 
before  returning  to  the  overlay  that  requested  a 
FACT  run. 

EGSDUCTiSR  The  main  Surface-Duct  Model  routine, 
called  from  FACTL3  as  SDUCT.  The  physical  and 
mathematical  basis  for  this  routine  and  its 
supporting  routines  is  documented  in  Reference  17. 

EGSDCOM:SR  A Block  Data  Subprogram  to  initialize 
common  area  variables  for  the  Surface-Duct  Model. 

EGSDINIT:  SR  The  Surface-Duct  input  parameter 
initialization  routine.  The  routine  computes  the 
equivalent  bilinear  (two  lines)  sound  speed 
profile  based  on  the  complete  sound  speed  profile. 
It  calls  subroutine  LSFIT  to  perform  a least 
squares  fit  to  the  in-layer  and  below-layer 
segments  of  the  profile 

EGLSFITtSR  This  routine  does  a least  squares  fit 
to  the  set  of  profile  depth  points  Z(l),  ....  Z(N) 
and  the  corresponding  functional  values 

F( I ) - (1  - C(l)/C(I)]2,  1-1,  ....  N 

where  C ( I ) are  the  profile  sound  speeds.  The 
function  was  chosen  because  the  Surface-Duct  Model 


asauroea..this  function  to  be  linear  for  the  in- 
layer and  below-layer  segments  of  the  profile. 

EGTLOUT;SR  This  routine  is  called  from  FACTL3 
to  display  the  TL  array  for  debug  purposes  only. 

EGSLOSS  :SR  FACT  rough  surface  reflection-loss 
routi.ce,  called  by  SDUCT. 

FACTLIB3  ;SR  This  library  of  routines  includes 
supporting  software  for  the  main  FACT  surface-duct 
routine,  SDUCT.  It  includes  LEAKY,  TRAPED,  SMODE, 
HISQ,  REFLK,  AIRY  and  FANG. 

Surface-Duct  Model  Support  Software 

FACT:  SR  A test  driver  program  for  the  FACT 
propagation-loss  model.  It  transfers  control  to 
the  FACT  overlays  ( FACTL1 , FACTL2 , and  FACTL3) , 
and  displays  the  results  numerically  and/or 
graphical ly . 

EGSDMAIN : SR  A test  driver  program  for  the 
Surface-Duct  Model.  For  comparison,  the  user  may 
select  either  the  old  model  (Clay)  or  the  new 
model  (SAI/Ryan).  Two  versions  of  the  program  can 
be  executed,  one  without  graphics  output,  but  with 
detailed  debug  output  capability  (EGSDMAIN);  and 
the  other  with  graphics  output  but  with  no  de- 
tailed debug  output  capability  ( EGSDPLOT) . 


EGSPINMjSR 


A test  driver  program  for  subroutine 

SDINIT. 

EGTLPLTsSR  A transmission-loss  graphics  output 
routine. 

3.2.3  Surface-Duct  Model  - Special  Considerations 

The  Surface-Duct  Model  uses  a number  of  complex 
arithmetic  relationships.  During  the  course  of  testing  the 
model  a number  of  inconsistencies  were  found  when  the  CDC 
6600  and  DEC-10  versions  of  the  model  were  compared  with  the 
TAPPS  Nova  800  implementation.  The  lack  of  agreement  was 
traced  to  a number  of  places  where  the  Data  General  com- 
piler did  not  generate  correct  code  for  FORTRAN  statements 
which  use  various  types  of  complex  arithmetic  expressions. 
Although  it  was  possible  to  program  around  these  problems  in 
all  cases,  care  must  be  taken  if  the  model  is  to  be  modified 
in  the  future  to  ensure  correct  results. 

3.3  AMBIENT  NOISE  SOFTWARE  MODIFICATIONS 

3.3.1  TAPPS  Ambient  Noise  Software  Modifications  \ 

"1—  ' " ' ■ | " ' ' « 

NSLIBrSR  (SHPSR)  Subroutine  SHPSR  was  modified 
to  handle  the  new  shipping-noise  file  (ANFILE), 
which  was  expanded  to  include  the  southern  hemi- 
sphere. 

ANFILE  This  file  was  modified  to  include  the 
Indian  Ocean.  Also,  provision  was  made  for  the 
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complete  southern  hemisphere  and  some  entries  in 
the  Mediterranean  Sea  were  found  to  be  in  error 
and  subsequently  modified. 

3.3.2  Ambient  Noise  Support  Software 

ANFCREAT ; SR  This  routine  creates  ANFILE.  It 
accepts  as  input  the  ASCII  file  ANFILEtSR,  which 
contains  omni-directional  shipping  levels  for  each 
five  degree  square. 

3.4  ROUGH-SURFACE  LOSS  SOFTWARE  MODIFICATIONS 

3.4.1  TAPPS  Rough-Surface  Loss  Software  Modifications 

FACTL2 : SR  The  variable  CSC2  was  added  to  compute 
the  ratio  of  sound  speed  at  the  surface  to  the 
sound  speed  at  C(K2).  CSC2  ara  WHF  (wave  height) 
are  passed  as  arguments  to  subroutines  CUSP  and 
INSTOR  for  rough  surface  calculations. 

FACTLIB2:SR  (CUSP)  CUSP  was  modified  to  compute, 
the  sin  squared  of  the  grazing  angle  at  the' 

» 

surface  and  to  call  SLOSS  for  rough-surface 
losses. 

FACTLIB2 : SR  (INSTOR)  INSTOR  was  modified  to  com- 
pute the  sin  squared  of  the  grazing  angle  at  the 
surface  and  to  call  SLOSS  for  rough-surface 
losses . 


I 

I 
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EGSDUCT:SR  Subroutine  SDUCT  was  modified  to  call 


SLOSS  for  the  rough-surface  loss  contribution  to 
the  surface-duct  propagation. 

EGSLOSS ; SR  Subroutine  SLOSS  was  created  to  com- 
pute the  rough-surface  loss  as  a function  of  wave 
height,  the  sin  of  the  surface  grazing  angle 
squared,  the  number  of  surface  reflections,  and 
the  frequency. 

3.4.2  Rough-Surface  Loss  Support  Software 

EGSLSMN : SR  A test  driver  program  for  subroutine 
SLOSS, 

3.5  MISCELLANEOUS  SOFTWARE  MODIFICATIONS 

NSLIB: SR  (SHIPMS)  The  variable  SCOR  was  limited 
to  be  at  least  lO-1^  to  prevent.  TAPPS  from  trying 
to  take  the  log  of  a negative  number.  This  is 
equivalent  to  limiting  the  shipping  noise  correc- 
tion level  to  be  at  least  -150  dB  on  each  beam. 

DBLIB;SR  (TLGET)  The  variable  NSAMP  (the  number 
of  depth  points  computed  in  PREOPT)  was  tested  to 
detect  the  case  where  PREOPT  does  not  compute  any 
depth  points. 

GETENV : SR  The  variable  SHPDEN  was  previously 
spelled  SHIPDEN  in  several  places. 


3-23 


INITsLD 


The  Assign  Console  File  INITrLD  did  not 
save  the  memory  locations  3542-3640,  which  con- 
tained initialization  for  the  variables  ALDZ, 
ADMIN,  EFUNC  and  AFUNC  (Common  area  ANAP) . These 
locations  were  added  to  the  save  command. 
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Section  4 
TEST  RUNS 


A number  of  test  cases  were  run  to  verify  the 
TAPPS  software  modifications.  Prior  to  full  integration 
within  TAPPS  each  of  the  major  subroutines  and  data  files 
were  tested  using  test  driver  programs: 


LSDBLP  - read  and  display  any  TAPPS  Bottom 

Loss  class  data  file  (e.g., 
INDABLP) 


LSBTM/LSDBLOS  - call  BTMLOS  to  compute  and 

display  bottom  loss  versus  grazing 
angle,  frequency  and  bottom  class. 
LSDBLOS  displays  results  in 
graphical  format. 

EGSDINM  - test  subroutine  SDINIT 


EGSDMAIN/ 

EGSDPLOT 


FACT 


EGSLSMN 


test  surface-duct  program  in  a 
stand-alone  mode  (subroutine 
SDUCT) , EGSDPLOT  displays  results 
in  graphical  format 

test  FACT  in  a stand-alone  mode 
(overlays  FACTL1,  FACTL2  and 
FACTL3 

test  rough-surface  loss  routine 
(SLOSS) 


The  source  (:SR)  and  loader  files  (:LD)  for  each 
of  the  above  test  programs  has  been  retained  on  the  TAPPS 
master  source  (fixed)  disk  for  additional  testing  if  re- 
quired . 


A number  of  full-scale  TAPPS  runs  were  also  run 
to  test  the  full  impact  of  the  model  changes. 
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Assign  Console  files  were  developed  for  many  of 
the  FACT  and  TAPPS  test  cases  and  have  been  retained  on  the 
TAPPS  master  load  (removable)  disk  for  additional  use  if 
needed.  The  files  are  named  FACTXXX:AC  and  TAPPSXX:AC. 

Figure  13  contains  a FACT  stand-alone  test  run, 
with  the  TAPPS  Version  1.2  modifications.  The  Assign 
Console  file  is  FACT331:AC. 

Figures  14a,  14b  and  14c  contain  TAPPS  test  runs 
for  TAPPS  Version  1.1,  with  the  BLUG  modifications,  and  for 
TAPPS  Version  1.2.  The  Assign  Console  file  is  TAPPS04:AC. 
These  test  runs  dramatically  illustrate  the  impact  of  the 
software  modifications  on  the  TAPPS  performance  estimates. 


\ 

I 
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T331  sac 


• P FACT 

PUMPER  OF  RANGE  INCREMENTS  (--.LE.250 ;-250 

pan  :»e  increment  i.nm)  «i 

GPmPH  LIMITS (MIN, MAX ) (DB  > -70.,  130. 

LIST, GRAPHICS  OUTPUT  (0-NO, 1-VES )-0, 1 
UAv'E  HEIGHT  (FT) -5 

LOU  FREQUENCV  BOTTOM  CLASS (1-3  OR  4-9)  -1 

HIGH  FREQUENCV  BOTTOM  CLASSC1-9)  -1 

SOURCE  DEPTH (FT) -60 

RECEIVER  DEPTH (FT  )-60 

NUMBER  OF  FREQUENCIES ( 1-5 ) -3 

FREQUENCY (HZ ) -50 

FREQUENCY (H2) -200 

FREQUENCY* HZ) -800  • 

INDEX  OF  LAYER -2 

NUMBER  OF  DEPTH-UELOCITY  PAIRS-5 

DEPTH^O 

VELOCITY -5046 

DEPTH-200 

UELOCITY-5050 

DEPTH- 1000 

VELOCITY -4960 

DEPTH-6000 

VELOCITY-4905 

DEPTH- 19000 

VELOCITY-5100 

DEBUG  MODE C 0-NO, 1 -YES >-0 

CALL  FACTL1 

SDINIT-IML, YS, YR, Z,C  t 4 50.0  ' 60, C 

0.0  50.0  60.0  200.0  1000.0  6000 ♦ 919008 ."6 

5046.0  5047. 0 5047.2  5050.0  4960.2  4906.4  5104.6 
LSFI T-KUT, ZD, N, 2, C s 0 200.001  4 

0.0  50.0  60.0  200.0 


Figure  13:  FACT  Test  Run  FACT  331: AC 
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Figure  13  (Continued):  FACT  Test  Run  FACT  331. ’AC 


Figure  13  (Continued):  FACT  Test  Run  FACT  331:AC 
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TAPPS  Test  Run  TAPPS04:AC  (Version 


TAPPS04 : AC  (Version 


(Continued):  TAPPS  Test  Run  TAPPS 04 : AC  (Version 


Figure  14b:  TAPPS  Test  Run  TAPPS04 : AC  (with  BLUG  Mods  Only) 
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Figure  14b  (ContinueaT:  TAPPS  Test  Run  TAPPS04:AC  (with  BLUG  Mods  Only) 


Figure  14b  (Continued):  TAPPS  Test  Run  TAPPS04:AC  (with  B LUG  Mods  Only) 
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Figure  14b  (Continued):,  TAPPS  Test  Run  TAPPS04:AC  (with  BLUG  Mods  Only) 
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Figure  14c  (Continued):^  TAPPS  Test  Run  TAPPS04:AC  (Version 
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